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This paper examines the response of electrolyte-saturated polymer gels, embedded with charged 
spherical inclusions, to a weak gradient of electrolyte concentration. These composites present 
a model system to study microscale electrokinetic transport processes, and a rigorous theoretical 
prediction of the bulk properties will benefit novel diagnostic applications. An electrokinetic model 
was presented in an earlier publication, and the response of homogeneous composites to a weak 
electric field was calculated. In this work, the influence of the inclusions on bulk ion fluxes and 
the strength of an electric field (or membrane diffusion potential) induced by the bulk electrolyte 
concentration gradient are computed. Effective ion diffusion coefficients are significantly altered by 
the inclusions, so — depending on the inclusion surface charge or (^-potential — asymmetric electrolytes 
can behave as symmetrical electrolytes, and vice versa. The theory also quantifies the strength of 
flow driven by concentration-gradient-induced perturbations to the equilibrium diffuse double layers. 
Similarly to diffusiophoresis, the flow may be either up or down the applied concentration gradient. 

I. INTRODUCTION 

Permeable membranes are the basis of molecular separation and sorting technologies, including ion exchange, 
gel-electrophoresis chromatography, and dialysisi^. Ion permeable membranes are also vital to batteries, fuel cells, 
electrochemical sensors, and biological cells^S. In this work, a theoretical model of electrolyte transport in membranes 
comprised of a continuous polymer gel with charged spherical inclusions is used to investigate the influence of the 
inclusions on bulk ion fluxes. Because several important characteristics of the mierostructure can, in principle, 
be carefully controlled, these materials present a model system for studying fundamental aspects of electrokinetic 
transport processes. An ability to tailor the mierostructure and accurately predict the bulk properties could also 
make these materials useful in novel diagnostic applications^. 

Electrolyte diffusion in the absence of electromigration is accompanied by a net flux of charge when anion and 
cation mobilities are not equal. Therefore, an electric field is necessary to compensate for the net diffusive flux of 
charge across a membrane. In biological cells, ion mobilities are regulated by transmembrane proteins (ion channels), 
and electrical signaling is achieved by controlling the relative ion mobilities'^. In diagnostic electrochemical cells, 
diffusion potentials are deleterious, and their influence is usually attenuated by a salt bridge containing a concentrated 
symmetrical electrolyte [e.g., KCl). When charged inclusions are immobilized in an electrically neutral polymer gel, 
the effective ion mobilities are altered, so even symmetrical electrolytes my behave as highly asymmetric electrolytes; 
this paper provides a first step toward quantifying this influence. 

When aqueous NaCl, for example, is brought into contact with negatively charged inclusions (e.i?., silica), the 
counter- ion (Na"*") is concentrated in the diffuse double layers. A gradient of electrolyte concentration perturbs the 
equilibrium state, inducing diffusion and electromigration to restore local equilibrium. At the surfaces of inclusions 
facing up the concentration gradient, the equilibrium gradient of electrostatic potential attracts counter- ions and 
repels co-ions, leading to inner and outer layers of perturbed charge density. The accompanying electric field acts 
on the underlying equilibrium charge, driving layers of backward and forward electroosmotic flow. In this work, 
non-linear perturbations are neglected because the perturbing 'forces' are weak compared to those of the underlying 
equilibrium state. Accordingly, the theory is limited to weak electrolyte concentration gradients. 

Perturbations enhance the net flux of counter- ions (Na"*"), because counter-ions are accumulated (depleted) at 
the surfaces facing up (down) the bulk concentration gradient. The resulting (tangential) concentration gradient 
increases the net (diffusive) counter-ion flux. Note that the accompanying perturbed charge density produces an 
electric field that induces an electromigrative flux of counter-ions down the bulk concentration gradient. Therefore, 
because Na"*" ions have a lower mobility than Cl~ , both contributions to the net flux counteract the tendency of Na+ 
to otherwise diffuse more slowly than Cl~. Overall, negatively charged inclusions increase the (effective) symmetry 
of the electrolyte. Since the degree of electrolyte symmetry bears directly on the membrane diffusion potential, the 
diffusion potential characterizes the physicochcmical state of the inclusion-electrolyte interface. 

The theoretical model adopted in this work was developed and applied in an earlier publication to examine ion fluxes 
and electroosmotic flow driven by an applied electric fieldS. In homogeneous membranes, diffusive and electromigrative 
transport are significantly infiuenced by the charge of the inclusions, but not by the hydrodynamic permeability of 
the polymer gel. On the other hand, electroosmotic flow is very sensitive to the gel permeability and reorganization 
of charge by diffusion and electromigration. These characteristics are expected to prevail with the application of a 




FIG. 1: Schematic of the microscale under consideration. Charged, impenetrable inclusions {solid circles) with radius a ^ 
10 nm-10 fim are embedded in a continuous polymer gel {solid filaments) saturated with aqueous electrolyte. Diffuse double 
layers {dash-dotted circles) with thickness k.^^ ~ 1-100 nm are perturbed by the application of an average electric field — (V-i/;), 
pressure gradient (Vp), or electrolyte concentration gradient (Vrij). The Brinkman screening length I ~ 1-10 nm that specifies 
the Darcy permeability £^ of the gel is often small compared to the inclusion radius. 



bulk concentration gradient, as examined in this paper. 

Note that there is a close connection of this work to earlier studies of diffusiophoresisiiii^'^, which address the 
motion of colloidal particles induced by solute concentration gradients. Fluid motion induced by solute concentration 
gradients at stationary charged interfaces is referred to as diffusioosmosi^, and this, perhaps, best characterizes the 
phenomena investigated here. Note that Wei and KehS^ derived semi-analytical predictions of the flow induced by 
solute concentration gradients in (ordered) fibrous porous media. Their calculations are restricted to low surface 
potentials, however, and their theory adopts a so-called cell model to handle fiber-fiber interactions. In the limit of 
vanishing fiber volume fraction, the fibers are intervened by pure solvent, whereas in this work the inclusions (spheres) 
are mediated by an uncharged Brinkman medium. 

The paper is organized as follows. The electrokinetic model and methodology are outlined in ^] Results demon- 
strating the response to a weak electrolyte concentration gradient, in the absence of a bulk electric field, are presented 
in ijlll Al These highlight the influence of negatively charged inclusions on the effective diffusivities of the electrolyte 
(NaCl) ions. Section IIII Al also examines the prevailing electroosmotic flow, which is either backward or forward, 
depending on the surface charge, gel permeability, and ionic strength. Superposition of two linearly independent 
problems is used in i jlll Bl to solve the problem with co-linear bulk electrolyte concentration and electrostatic po- 
tential gradients. With the constraint of zero bulk current density, the ratio of the electric field strength to the 
concentration gradient is established. A brief sumnrary follows in illVI 



II. THEORY 



The microstructure of the composites addressed in this work is depicted schematically in figure ^ The continuous 
phase is a porous medium comprised of an electrically neutral, electrolyte-saturated polymer gel. Polyacrylamide 
gels are routinely used for the electrophoretic separation of DNA fragments in aqueous media. The porosity is 
controlled by adjusting the concentrations and ratio of the monomer (acrylamide) and cross-linker. In this work, the 
hydrodynamic permeability is characterized by the Darcy permeability £^ (square of the Brinkman screening length), 
which reflects the hydrodynamic size and concentration Ug of the polymer segments. In turn, these reflect the 
degree of cross-linking and the affinity of the polymer for the solvent. 

Embedded in the polymer are randomly dispersed spherical inclusions. In model systems, these inclusions are 
envisioned to be monodisperse silica beads or polymer latices, with radii in the range a = 10 nm-10 fim. The 
inclusions bear a surface charge when dispersed in aqueous media, and the surface charge density may vary with the 
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bulk ionic strength and pH of the electrolyte. In this work, the charge is to be inferred from the bulk ionic strength 
and electrostatic surface potential, ^. Because the inclusions are impenetrable with zero surface capacitance and 
conductivity, no- flux and no-slip boundary conditions apply at their surfaces. 

Note that the mobile ions whose charge is opposite to the surface-bound immobile charge are referred to as counter- 
ions, with the other species referred to as co-ions. For simplicity, the counter-charges (dissociated counter-ions) are 
assumed indistinguishable from the electrolyte counter-ions. Surrounding each inclusion is a diffuse layer of mobile 
charge, with Debye thickness k^^ and excess of countcr-ions; as described below, the equilibrium double-layer structure 
is calculated from the well-known Poisson-Boltzmann equation. 



A. Microscale model 



The microscale transport equations and boundary conditions are based on the standard electrokinetic mode&^ with 
a body force to model the frictional resistance of the polymes^i^. This coupling has been widely used to interpret the 
electrophoretic mobility and other characteristics of 'soft' colloidal particles and their dispersionsiSiiiiSiSiiSi. The 
model comprises the (non-linear) Poisson-Boltzmann equation 

N 

eoe^V'^ip ^ -^UjZjC, (1) 

where and are the permittivity of a vacuum and dielectric constant of the electrolyte; Uj are the concentrations of 
the jth mobile ions with valences Zj ; and ip and e are the electrostatic potential and elementary charge. Ion transport 
is governed by the Nernst-Plank relationship 

67r7^aj(M - Vj) - ZjcVip - fcTV Inuj =0 {j ^ 1, TV), (2) 

where Oj are the Stokes radii of the ions, obtained from limiting conductances or diffusivities; rj is the electrolyte 
viscosity; u and Vj are the fluid and ion velocities; and kT is the thermal energy. Ion diffusion coefhcients, which will 
be adopted throughout, are 

Dj = kT/{6TTT]aj). (3) 
As usual, the double-layer thickness (Debye length) 



K-^ = v/A:Te,eo/(2/e2) (4) 



emerges from Eqns. ^ and (|2Jl where 



N 

I ={1/2)Y, z]n^ (5) 
i=i 

is the bulk (average) ionic strength, with the bulk ion concentrations. 
Ion conservation demands 

dn,/dt = Q=-V-{j^) (j = l,...,iV), (6) 

where t is the time, with the ion fluxes jj ~ n-jVj obtained from Eqn. jSJ. Similarly, momentum and mass conservation 
require 

N 

Psdu/dt = = TyV^M - Vp - (?//^^)m - ^ UjZjeVi) (7) 

and 

V • M = 0, (8) 

where ps and u are the electrolyte density and velocity, andp is the pressure. Note that —(r]/£'^)u is the hydrodynamic 
drag force exerted by the polymer on the electrolyte. The Darcy permeability (square of the Brinkman screening 
length) of the gel may be expressed as 

= l/[n,{r)67ra,Fs] = 2al/[9Mr)Fs{<j),)], (9) 
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where ns{r) is the concentration of Stokes resistance centers, with Os and Fg the Stokes radius and drag coefficient 
of the polymer segments. The hydrodynamic volume fraction 0^ = ns{4:/i)'!Tal of the segments in swollen polymer 
gels is often very low, so Fs « 1. In this work, the Brinkman screening length is adjusted according to Eqn. Q by 
varying the (uniform) polymer segment density with Stokes radius = 1 A. The drag coefficient is obtained from 
a correlation for random fixed beds of spheresii. For the purposes of this paper, however, only the reported values of 
£ are relevantf*^. 

Either the equilibrium surface potential C or surface charge density a may be specified. Because the surface (r = a) 
is assumed impenetrable with zero capacitance and conductivity, the surface charge is constant, permitting no-flux 
boundary conditions for each (mobile) ion species. As usual, the no-slip boundary condition applies. 

In the far field, neglect of particle interactions requires 

ip —E ■ r as r ^ oo, (10) 



rij n°° + Bj • r as r oo, (11) 

and 

u ^ U as r ^ oo, (12) 

where E, Bj and U are, respectively, a constant electric field, constant species concentration gradients, and constant 
far-field velocity. 
With 'forcing' 

X = Xe,, (13) 

where X E {E,Bj,U}, linearized perturbations —E-r + tl'', Bj ■ r + n'j and U + u' to the equilibrium state (-0°, n^j 
and tt = 0) are symmetric about the z-axis {6 = 0) of a spherical polar coordinate system. Primed quantities have 
the form 

il)' = i}{r)X ■ Br (14) 



n'j ~ hj{r)X ■ Er (15) 

and 

u' = V X V X h{r)X 

= -2{hr/r){X ■er)er-{hrr + hr/r){X ■eg)ee. (16) 

The radially decaying functions il>{r), nj{r) and h{r) are calculated numerically in this work. However, as discussed 
below, and at length elsewhere^, only the scalar coefficients characterizing their far-field decays (see Eqns. ((THt - (pn|) 
below) are necessary to derive bulk properties of the composite. Therefore, the practical purpose of the numerical 
procedure is to determine these so-called asymptotic coefficients. 

Solutions of Eqns. (J^), ® > O a-nd l|SJ), to linear order in perturbations to the equilibrium state, with E, Bj and U 
set to arbitrary values, can be computed if 

N 

Y^ZjB.^O (17) 

to ensure an electrically neutral far-field. However, when N species are assembled into M clcctroneutral groups {e.g., 
electrolytes or neutral tracers), each with far-field gradient Bk {k = 1,...,M), it is expedient to compute solutions 
with only one non-zero value of E, Bk or U. Then, arbitrary solutions can be constructed by linear superpositioni^. 

An index k' is required to identify the (electroneutral) group to which the jth species under consideration is 
assigned. Careful consideration of the electrolyte composition and ion valences is required to ensure consistency. For 
z-z electrolytes it is convenient to set Bj = Bk, whereas for a single 1-2 electrolyte {e.g., CaCl2) is it satisfactory to 
set Bj = Bi;' /\zj\. For the relatively simple situations considered in this work, the (single) electroneutral group is 
NaCl, so M = 1 with k = k' = 1, and j = 1 and 2 for Na+ and Cl^, respectively. 
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The perturbations satisfy a linear set of coupled ordinary differential equations with far-field boundary conditions^ 

-^{X ■ er)D^ Ir^ as r oo, (18) 



n^. ^ (X • e.r)Cf jr"^ as r ^ oo, (19) 

and 

m' ^ -2(C^/r3)(X • er)er + (C^ lr^){X ■ eg)eg as r ^ cx). (20) 

The equations for the non-linear equilibrium state and the linearized perturbations are solved using the numerical 
methodology of Hill, Saville and Russeliii, which was developed for the electrophoretic mobility of polymer-coated 
colloids and the bulk properties of their dilute dispersions. 

Note that C^*" is the asymptotic coefficient for the perturbed concentration of the jth species induced by the 
fcth concentration gradient Bj^, whereas C^*" (without a subscript) denotes the asymptotic coefficient for the flow 
induced by Bk- For neutral species, the concentration disturbance produced by a single impenetrable sphere yields 
C^*"' — {l/2)a^, otherwise C^'' = (fc 7^ k'). Clearly, the asymptotic coefficients for charged species, whose 
concentration perturbations are influenced by electromigration. arc not the same as for neutral species; for ions, 
however, cf"' (l/2)a3 as |C| ^ 0. 

With co-linear forcing and bulk electroncutrality, linear superposition gives far-field decays 

M 

ip' -> {l/r^)[ED'^ + J2BkD^'' + UD^]{e,^ ■ e,.) as r ^ 00, (21) 

k=l 
M 

n; {l/r'')[ECf + BkCf" + UC^]{e, • e,) as r ^ ^, (22) 
fc=i 

and 

M 

u' -{2/r^)[EC^ + Y,BkC^'' +UC'^]{e^-er)er 

k=l 

M 

+ {l/r^)[EC'' + BkC^' + UC^]{e, ■ eg)eg as r ^ 00. (23) 

fc=i 

These relationships are the key results from which all macroscale quantities {e.g., bulk ion fluxes and electroosmotic 
flow) are derived for small inclusion volume fractional. Situations with only one forcing variable X G S^, ?7} are 
referred to as the (E), (B) and (U) (microscale) problems. Algebraic or differential relationships between the averaged 
fields can be applied to ensure, for example, zero average current density (see IIII B|) . 

B. Macroscale equations 

Expressions relating the asymptotic coefficients , and to the averaged ion fluxes and momentum conser- 
vation equation can be derived using procedures similar to those applied by Saville^Sand O'Brieni^for the conductivity 
of colloidal dispersions. As demonstrated in an earlier paper—, ion and fluid momentum fluxes can be averaged over 
a representative elementary control volume, and the volume integrals enumerated from knowledge of the asymptotic 
coefficients. Similarly to Maxwell's well-known analysis of conduction, particle interactions are neglected, so the re- 
sults arc limited to small volume fractions (f)[\ + {Ka)~^Y ^ 1? where (jj — ri,(4/3)7ra'^ is the inclusion volume fraction, 
with n the inclusion number density. 

When all average fluxes are in the z-dircction, mass and momentum conservation require constant average velocity 
(■u), giving^ 

M 

(Vp) = -(r//£2)(M) - ^{i/a^){T^ie)[EC^ + ^ BkC^" + UC% (24) 

k=l 
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where (Vp) is the average pressure gradient (set to zero in this work). Similarly, the (steady) average species 
conservation equations V • {j j) = require constant average fluxes^ 

n *^ 

+(b{3/a'')zje^n°°[ED^ + ^ BkD^'^ + UD^] 

k=l 
AI 

+^{3/a^)Dj[ECf + BkCf" + JJCf]. (25) 

k=l 

Note that V- (Vi/^) = in an electrically neutral composite with uniform dielectric permittivity, so the average electric 
field (V'i/') is also constant. 

The averages can be expanded as power series in the inclusion volume fraction e.g., (u) Uq + 4>Ui + 0{(j)^). 
Therefore, since the microscale equations (asymptotic coefficients) are accurate to 0((j)), the notation is condensed 
by writing, for example, (it) = U, where it is understood that U = Uq + 4>Ui + 0{(jP'). Clearly, E, Bj and U in 
Eqns. (|24|l and H25|l need only include the 0(1) contribution to their respective average value, e.g., U — > Uq. The 
following notation is adopted for the other averaged quantities: Jj = (jj), P = (Vp), Bj = (Vrij), E = — (VV'). 

With one electrolyte (Af = 1) and, recall, bulk electroncutrality, there are iV + 4 independent variables 
{E,U ,P,Bk {k = 1), Jj (j = l,...,iV)) with iV + 1 independent equations (Eqns. ij^H) and (ESJ)- Clearly, three 
independent variables must be specified for a unique solution. 

For clarity, the results presented below involve a 1-1 electrolyte (NaCl). It is important to note that, because the 
equations are linear, solutions for any combination of non-zero forcing variables may be constructed. For example, 
^III BI establishes the electric field strength required to maintain a constant electrolyte flux — driven by a bulk con- 
centration gradient across a membrane — with zero bulk current density. Solutions of the (E) problem are available 
elsewhere^, whereas all solutions of the (B) problem presented below are new. 

III. RESULTS 
A. Concentration gradient alone 

When an average concentration gradient is applied in the absence of an average pressure gradient and electric field, 
the ion fluxes (Eqn. (|25|) ') are 

n ^ 

J, = -D,B, + 0(3/a3)z,e^nf ^ B^D^^ 

k=l 

M M 

+^{3/a^)D, J2 BkCf" + (/.(3/a3)nf ^ B^C^" + 0(0^). (26) 

k=l k=l 

The first term on the right-hand side is the diffusive flux in the absence of inclusions, and the second, third and fourth 
terms, respectively, are corrections due to microscale electromigration, diffusion and convection due to the inclusions. 

For example, with the composite bridging two reservoirs, one containing KCl with concentrations at position 
z = 0, and the other containing NaCl with concentrations n°°2 at z = L, it is necessary to designate M = 2 bulk 
electrolyte concentration gradients for the = 3 species: i?i = —n^-^/L for KCl, and B2 = n^2l^ foi' NaCl. 
Therefore, the bulk concentration gradients for K+, Na+, and Cl~ are, respectively, B\, B2 and Bi + B2. 

For the simple case involving the concentration gradient of a single z-z electrolyte {Bj = Bk', iV = 2 and M = 1), 
Eqn. H26|) may be written 

J, =-i?,B,(l + <^Af'), (27) 



where 



-{i/a^)^y^D^-' - i3/a^)Cf'" - (3/a^)-^C^'^' . (28) 
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This motivates the introduction of an effective diffusivity 

D* = - J,/B, =D,{1+ 0Af ) , (29) 

where Aj^ = A^*"' is termed the (dimensionless) effective diffusivity increment. 

Note that electromigration from concentration-gradient-induced electrical polarization increases (decreases) the 
bulk counter-ion (co-ion) flux, leading to a (bulk) current density 

N 

N 

= -Bk'J2z,eD,{l + cl)Af'''). (31) 

An average concentration gradient also induces an 0{4>) average velocity that reflects the permeability i'^ of the 
polymer gel. ft is therefore expedient to examine the ratio of the average fluid velocity to the product of the particle 
volume fraction and the strength of the concentration gradient 

U/iBk'(f>)^-3C''^'/a\ (32) 

This quantity, termed the incremental pore mobility, is similar to the incremental pore mobility that characterizes 
electric-field-induced flow-. It is also bears similarity to the diffusiophoretic mobility of free colloidal particlesiS. 

Note that Eqns. (|31|l and H32|l apply with zero average electric field. This implies that an electric field is impressed 
to counteract the field that would otherwise develop spontaneously to achieve zero bulk electrical current. The pore 
mobility with zero electrical current, for example, may be obtained by eliminating E from Eqns. (|24|l and (|25|) . with 
P = I = 0, as demonstrated in Wl Bl For the simple case involving a single z-z electrolyte {Bj = B^', N — 2 and 
A/ = 1), the pore mobility is 

U/{Bk'(t>) - -{3/a^)[C^''' +C^J2 + 0{(f>), (33) 

where 

=Y,{z,ef^nf (34) 
i=i 

is the electrolyte conductivity. 

Asymptotic coefficients are provided in table for a representative polymer gel with Brinkman screening length 
£ « 0.951 nm and inclusion radius a ~ 100 nm; the C-potentials and (three) ionic strengths span experimentally 
accessible ranges. Here, the sign of the electrostatic and concentration dipole moments (D^'-' and Cf*'' , respectively) 
and of the far-field flow (as indicated by (7^*=' ) are not straightforward to interpret without examining the detailed 
structure of the perturbed double layer. Representative (dimensionless) functions fij and ■0(r) are shown in figure 13 
(see Eqns. (I18|) and H19() ). Recall, these characterize the radial variation of the perturbed ion densities and electrostatic 
potential. Figure |2| shows the corresponding radial fluid velocity, as represented by —2hr/r (see Eqn. (EHJ)- 

From table backward concentration-gradicnt-induced flow {U < or C^''' > 0) prevails at low to moderate 
values of na, with forward flow {U > 0) evident at higher values of kg when \(\ is low. Recall, in the absence of 
surface charge (^ = 0), the concentration dipole moments C^^*"' = (1/2)0"^ reflect the (concentration) polarization 
necessary for electrolyte to diffuse past impenetrable inclusions. With surface charge, however, counter-ions (co-ions) 
at the forward facing surfaces migrate radially inward (outward) under the influence of the equilibrium electrostatic 
potential. The resulting inner and outer layers of perturbed charge density are evident in figure [5J these clearly have 
positive and negative dipole moments, respectively. Symmetry and electroneutrality considerations indicate that the 
sign of the net electrostatic dipole moment {i.e., the far-field decay of the electrostatic potential) reflects the moment 
of the outer layer of perturbed charge. Accordingly, the electrostatic dipole moments D^fc' in tableware negative for 
all values of kg and 

It is tempting to attribute the direction of the far-field flow to the sign of the electrostatic dipole moment, whose 
accompanying electric field is expected to drive forward flow when acting on the equilibrium charge density. However, 
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scaled radial distance, K(r-a) scaled radial distance, K(r-a) 



FIG. 2: The structure of the perturbed double layer — under the influence of a bulk electrolyte concentration gradient — as 
revealed by the (scaled) radial electrostatic potential, i/j(r)2/e/(fcTa) {dashed lines), and ion concentration, h\[r)/a (Na"*") and 
n2{r)/a (Cl^) [solid lines), perturbations (see Eqns. I|14^ - H16|l ') as a function of the (scaled) radial position K(r — a): aqueous 
NaCl at T = 25°C; a = 100 nm; I ^ 3.11 nm; (ko, Ce/(A:r)) = (1, -2) (left panel); {Ka,(e/{kT)) = (100, -2) (right panel). 

as indicated in tabled the far-field flow is most often backward at low to moderate values of kg. Furthermore, the 
direction of the far- field flow is independent of the sign of the net dipole moment. 

In situations where the far-field flow is backward, the inner perturbed electric field acts on the equilibrium charge 
density, and viscous stresses from the prevailing electroosmotic flow propagate to drive a backward outer flow. In 
this case, all streamlines are open, as indicated by the (negative) radial velocity profile in the left panel of figure |31 
However, as seen in the right panel, situations with a forward far-field flow also have a thin layer of (relatively slow) 
backward flow at the inclusion surface. In this case, closed streamlines prevail within the (thin) equilibrium double 
layer. Evidently, electrical forces due to the outer perturbed electric field acting on the equilibrium charge density 
dominate viscous stresses arising from the innermost layer of backward flow. 

While the qualitative features of these flows bear similarities to those underlying the diffusiophoretic motion of 
charged colloidal spheresiS, the differences are, perhaps, sufficient to discourage direct comparisons. For example, 
as seen below, the direction of the far-field flow in the composites addressed here also depends on the Brinkman 
permeability of the intervening polymer. Furthermore, the diffusiophoretic mobility is also influenced by an electric 
field (to maintain zero current density) whose sign depends on electrolyte asymmetry. 

For a composite with a relatively large Brinkman screening length i w 9.95 nm, the conccntration-gradient-induccd 
pore velocities shown in figure 0] are low. As expected, the incremental pore mobility increases with the surface 
charge (or ^-potential) and Darcy permeability. The mobility is shown in figure as a function of the Brinkman 
screening length for various ^-potentials, with an ionic strength / « 0.0925 mol 1~^ yielding ko = 100. The direction 
of the concentration-gradient-induced flow is consistent with the discussion of table U Increasing the gel permeability 
strengthens the (forward flowing) electroosmotic flow in the outer region of the perturbed double layers and, hence, 
increases the bulk convcctive flow. 

Contributions to the incremental fluxes are provided in table ^ for the composite whose asymptotic coefficients 
are listed in table ^ The first term on the right-hand side of Eqn. 128|) (also columns 2 and 3 of table is due 
to electromigration arising from concentration-gradient-induced electrical polarization. Polarization of (negatively) 
charged inclusions evidently enhances the effective diffusivity of the counter- ion (Na+) and attenuates that of the 
co-ion (Cl~). The second term in Eqn. (|28|) (column 4 of the table) is due to diffusion, and arises from perturbations 
to the ion concentration gradients. As expected, this attenuates the diffusive flux when the ^-potential is moderate, 
and approaches the Maxwell value A^^' = —3/2 as |C| 0. The last term in Eqn. H28|) (also column 5) is due 
to electroosmotic flow (convection). This diminishes the fluxes of both ions, but its contribution is relatively small. 
Summing all three contributions (last two columns) shows that the average diffusive flux of the co-ion (Cl^) is 
attenuated by the inclusions, most significantly at low ionic strength when electrostatic screening of the surface 
charge is weak. Similarly, the flux of the counter-ion (Na"*") is enhanced because of electromigration. Again, when 
electrical forces are weak, the inclusions simply hinder diffusion and —3/2 as |C| — > 0. 



TABLE L Asymptotic coefficients (see Eqn. 12611 1 for the bulk diffusion of NaCl in a Brinkman medium witli charged spherical inclusions: a = 100 nm; (. ~ 0.951 nm; 
T = 25°C, Di « 1.33 X lO'V^s-i (Na+); Da ~ 2.03 x lO-'^m^s"^ (Cr); «* = eseo{kT / ef / {-qa) « 5.15 x 10"^ m s'^ 
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-2 


-6.04 X 10"^ 


+3.98 X IQ-i 


9.99 X 10-^ 


-8. .4 X 10° 


-2.52 X lOi 


-4 


-1.29 X 10" 


+1.42 X IQ-i 


3.06 X 10-" 


-2.56 X lOi 


-5.37 X lOi 


-6 


-1.88 X 10" 


-1.20 X IQ-i 


4.12 X 10-1 


-3.44 X lOi 


-7.84 X lOi 


-8 


-2.23 X 10" 


-2.81 X IQ-i 


4.26 X 10-1 


-3.56 X lOi 


-9.29 X lOi 


= 100 


/ = 9.25 X 10"^ mol 










-1 


-3.08 X 10"^ 


4.96 X 10-' 


-3.01 X 10-** 


+2.52 X 10-' 


-1.28 X 10-2 


-2 


-6.85 X 10-2 


4.85 X 10-1 


-7.71 X 10-1 


+6.44 X 10-1 


-2.85 X 10-2 


-4 


-1.97 X 10-1 


4.28 X 10-1 


-6.87 X 10-* 


+5.74 X 10-1 


-8.19 X 10-2 


-6 


-4.49 X 10-^ 


3.04 X 10-1 


-3.05 X 10-^ 


+2.55 X 10-2 


-1.87 X 10-1 


-8 


-8.24 X 10-1 


1.18 X 10-1 


+3.15 X 10-1 


-2.63 X 10-1 


-3.43 X 10-1 



TABLE II: Scaled asymptotic coefficients and microscale incremental contributions (see text) to average diffusion of NaCl in a Brinkman medium with charged 
spherical inclusions: a = 100 nm; i fa 0.951 nm; T = 25°C; Di ^ 1.33 x 10"^m^s"^ (Na+); D2 « 2.03 x lO'^m^s"^ (CI"). 



(e/ikT) 


A"X (Na+) (= -A^f ) 


Km (J = 1.2) 


(Na+) 


A"r (cr) 


a""' (Na+) 


A^' (CI") 


na — 1 


/ = 9.25 X 10"'' mol 












-1 


6.25 X 10" 


-1.15 X 10° 


-8.76 X lO"*" 


-5.74 X lO"*" 


5.10 X 10° 


-7.39 X 10" 


-2 


1.22 X 10^ 


-1.71 X 10"' 


-3.02 X 10"^ 


-1.98 X 10"" 


1.21 X 10' 


-1.24 X 10' 


-4 


2.23 X 10^ 


+2.60 X 10° 


-6.61 X 10"'' 


-4.34 X 10"" 


2.49 X 10' 


-1.97 X 10' 


-6 


2.86 X 10^ 


+4.79 X 10° 


-6.61 X 10"^ 


-4.33 X 10"" 


3.34 X 10' 


-2.39 X 10' 


-8 


3.17 X 10^ 


+5.89 X 10° 


-5.66 X 10"^ 


-3.71 X 10"" 


3.75 X 10' 


-2.58 X 10' 


Ka = 10 


/ = 9.25 X 10"" mol 1"' 












-1 


4.38 X 10"' 


-1.42 X 10° 


-1.47 X 10"" 


-9.65 X 10"'' 


-9.84 X 10"' 


-1.86 X 10° 


-2 


9.07 X 10"^ 


-1.19 X 10° 


-5.79 X 10"" 


-3.80 X 10"" 


-2.88 X 10"' 


-2.10 X 10° 


-4 


1.94 X 10° 


-4.25 X 10"' 


-1.77 X 10"* 


-1.16 X 10"* 


+1.51 X 10° 


-2.36 X 10° 


-6 


2.82 X 10° 


+3.61 X 10"' 


-2.39 X 10"" 


-1.57 X 10"* 


+3.19 X 10° 


-2.46 X 10° 


-8 


3.35 X 10° 


+8.42 X 10"' 


-2.47 X 10"* 


-1.62 X 10"* 


+4.19 X 10° 


-2.50 X 10° 


Ka = 100 


/ = 9.25 X 10"^ mol 1"' 












-1 


4.63 X 10"^ 


-1.49 X 10° 


+1.75 X 10"* 


+1.14 X 10"* 


-1.44 X 10° 


-1.54 X 10° 


-2 


1.03 X 10"' 


-1.45 X 10° 


+4.47 X 10"* 


+2.93 X 10"* 


-1.35 X 10° 


-1.56 X 10° 


-4 


2.95 X 10"' 


-1.28 X 10° 


+3.98 X 10"* 


+2.61 X 10"* 


-9.88 X 10"' 


-1.58 X 10° 


-6 


6.74 X 10"' 


-9.13 X 10"' 


+1.77 X 10"" 


+1.16 X 10"" 


-2.39 X 10"' 


-1.59 X 10° 


-8 


1.24 X 10° 


-3.53 X 10"' 


-1.82 X 10"* 


-1.20 X 10"* 


+8.83 X 10"' 


-1.59 X 10° 
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scaled radial distance, K(r-a) scaled radial distance, K(r-a) 



FIG. 3; The structure of the perturbed double layer — under the influence of a bulk electrolyte concentration gradient — as 
revealed by the (dimensionless) radial velocity —2{hr/r)2I/{u*a) {u = eseo(kT/ef/{rta) « 5.15 x 10"^ m s"^) (see Eqn. (EgJ) 
as a function of the (scaled) radial position nir — a). Parameters are the same as in figure |21 




FIG. 4: The incremental pore mobility U/{By(j)) — —'iC^''' /a"^ as a function of the (scaled) ^-potential C^e/{kT) for (scaled) 
reciprocal double-layer thicknesses Ka = 1, 2 and 4: aqueous NaCl at T = 25°C; a = 100 nm; I « 9.95 nm. An average 
concentration gradient is applied in the absence of average pressure and electrostatic potential gradients. 

EfFcctive diffusivity increments arc shown in figure El as a function of the ^-potential for various values of na. 
Clearly, increasing the surface charge increases the effective diffusivity of the counter-ion (left panel) and decreases 
the diffusivity of the co-ion (right panel). Note that the calculations are insensitive to the gel permeability, so these 
results are applicable to a variety of composites with negatively charged inclusions and NaCl electrolyte. In general, 
all bulk properties that reflect concentration and electrical polarization are independent of the Darcy permeability. 

B. Bulk diffusion and electromigration 

Having examined the application of an electrolyte concentration gradient in the absence of an electric field, let us now 
consider the electric field strength and concentration gradient that together yield zero current density. This ensures 
that transport of ions from one reservoir to another (with different electrolyte concentration) maintains electrical 



12 




FIG. 5: The incremental pore mobility U/{Bki(f>) — — SC^^'/a^ as a function of the Brinkman screening length t for various 
(scaled) (^-potentials —C,e/{kT) = 1, 2, 4 {solid lines) 6 and 8 [dashed lines): aqueous NaCl at T = 25°C; a = 100 nm; na = 100; 
/ ?3 0.0925 mol r\ An average concentration gradient is applied in the absence of average pressure and electrostatic potential 
gradients. 




FIG. 6: The (scaled) effective diffusivity increments lS.f — {D*/Dj — l)/4> for Na^ [left panel) and CP {right panel) as a 
function of the (scaled) (^-potential (e/{kT) with various (scaled) reciprocal double-layer thicknesses na = 1, 2, 10, 100 and 
1000 (aqueous NaCl at T = 25° C with a = 100 nm). Note that results are insensitive to the permeability £^ of the polymer gel. 



neutrality across a membrane. As expected from the previous section, an electric field is necessary to compensate for 
the lower (higher) diffusive flux of the less (more) mobile ion. 

The analysis below is limited to situations where the bulk electrolyte concentration gradient is weak, permitting 
macroscale variations in the bulk equilibrium electrolyte concentration to be neglected. Complications that arise from 
macroscale variations in ion density, which are necessary to satisfy (bulk) continuity and electroneutrality constraints, 
will be addressed in future work; these arc discussed briefly in i llVI 

For the simple case involving the concentration gradient of a single z-z electrolyte {Bj ~ Bk' , N ~ 2 and AI = 1), 



13 



setting the current density (Eqn. (|30() ') to zero gives 

N=2 

EK°" = Bk' ZjeDj{l + <^A^) + O(02), (35) 

where 

« -(3/a^)[D^- ^ + Cf ' + i?^ + Cf ^^'-!.^ ^ ] + 0(0) (36) 

is termed the (dimensionless) electric-field increment. Note that and C^''' are independent of j, so is 
independent of the sum in Eqn. (|35|l . 

Similarly to the concentration-gradient-induced current density (Eqn. 1201, the 0(1) average electric field is zero 
for z-z electrolytes whose ions have equal mobilities. It is important to note that the 0(0) contribution is not zero, 
however, even for perfectly symmetrical electrolytes. This follows from the influence of charged inclusions on the 
effective ion diffusivities, and emerges from Eqn. (|36|l through the term involving the gradient-induced electrostatic 
dipole strength D^''' . Therefore, for a z-z electrolyte whose ions have equal diffusivities Dj, Eqn. (|35|l simplifies to 

E = -0(3/a3)i?^^-'Bfe, + O(02), (37) 

showing that the average electric field reflects a sum of concentration-gradicnt-induced electrostatic dipole moments. 
In an homogeneous membrane, the membrane diffusion potential A-0 = -0(2 = L) — "0(0) = —EL ~ (j){?>/a?)D^i'' B^iL 
is simply proportional to the concentration differential Bk'L. 

The electric-field increment is shown in figure as a function of the ^-potential for various bulk ionic strengths 
yielding Ka in the range 1-10"^. Again, because the increment reflects electrical and concentration polarization, the 
calculations are insensitive to flow and, hence, Darcy permeability. Therefore, the results are applicable to a variety 
of composites with negatively charged inclusions and NaCl electrolyte. Note that the increment is large at low ionic 
strength (small na) when the C-potential is high, so the particle contribution may be comparable to or greater than 
that of the electrolyte alone. Further, when A^ < 0, there exists an inclusion volume fraction 

0* = -1/A^ (38) 

that yields zero electric field strength with any (weak) bulk electrolyte concentration gradient. Clearly, 0* must be 
less than one to be consistent with assumptions underlying the theory. In turn, this limits the accuracy of such a 
prediction to situations where — A^ ~ 0~^ (0 ^ 1), which is achieved at low ionic strength {e.g., ko = 1 in figured. 

IV. SUMMARY 

This paper demonstrates that a weak electrolyte concentration gradient applied across a permeable membrane 
hosting immobile charged inclusions induces bulk clcctroosmotic (or "diffusioosmotic" ) fiow. The direction of this 
fiow may be up or down the concentration gradient, depending on the surface charge, bulk ionic strength, and 
permeability of the polymer gel. The inclusions also influence the ion fluxes or, equivalently, the effective ion diffusion 
coefficients. In short, the net counter-ion flux is enhanced by electromigration within the diffuse double layers. 

Steady or quasi-steady transport across a membrane often takes place with zero current density, so, in practice, an 
external electric field — or an electric field from a macroscale redistribution of space charge — is necessary to maintain 
electrical neutrality. This work established the proportionality between this electric field and the imposed concentra- 
tion gradient. It is tempting to attribute the strength of the electric field to electrolyte asymmetry alone, concluding, 
for example, that a perfectly symmetrical electrolyte will not produce a membrane diffusion potential. However, it 
was demonstrated that charged inclusions break the (effective) symmetry, so even a perfectly symmetrical electrolyte 
will yield a membrane diffusion potential. For NaCl (a moderately asymmetric z-z electrolyte) in the presence of 
negatively charged inclusions, the electric-field increment is negative, so the diffusion potential of the composite is 
lower than in the absence of inclusions. 

It should be emphasized that the present theory is limited to homogeneous membranes. Therefore, with a macro- 
scopic length scale L {e.g., membrane thickness), the differential concentrations An°° = BjL should be small relative 
to the bulk concentrations n°°. Even in the absence of inclusions, the bulk electrolyte fluxes are not straightforward 
to calculate without adopting one of the two following approximations. One consequence of electrical neutrality is 
that the electric field {e.g., with unidirectional transport) must be uniform. Similarly, ion conservation at steady 
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FIG. 7: The (scaled) electric-field increment (see Eqns. l).'-{5^ - H36|l 'l as a function of t lie (scaled) (^-potential —C,e/{kT) for 
various (scaled) reciprocal double-layer thicknesses na = 1, 2, 10, 100 and 1000: aqueous NaCl at T = 25°C; a = 100 nm. Note 
that the results are insensitive to the permeability of the polymer gel. 



state demands constant ion fluxes. However, solving the ion transport equations under these conditions leads to ion 
concentration fields that do not yield an electrically neutral bulk. This is the scenario that underlies the Goldnian- 
Hodgkin-Katz (GHK) equation, which is well known in membrane biologjf^^. An alternative approach constrains the 
ion conceirtration fields to ensure electrical neutrality. Then, however, the (bulk) ion fiuxes are not uniform and, 
hence, the ion concentrations cannot be steady. This scenario underlies the Henderson equation, which is well known 
in electrochemistry and electrochemical engineering^^. It is interesting to note that both theories yield the same 
membrane diffusion potential A?/' = ^{z = L) — "0(0) = —EL when the equations are linearized for small electrostatic 
potentials ip < kT/e. Given the difficulties in solving coupled electromigration and diffusion (electro-diffusion) in 
the absence of inclusions, the problem with charged inclusions seems intractable at present. This work will hopefully 
stimulate further work in this area. 
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